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Abstract. The leading-order accurate description of tt production, as usually employed in standard Monte Carlo 
event generators, gives no rise to the generation of a forward-backward asymmetry, Aps- Yet, non-negligible 
- differential as well as inclusive - asymmetries may be produced if coherent parton showering is used in the 
hadroproduction of top quark pairs. In this contribution we summarize the outcome of our study [ ! ] of this 
effect. We present a short comparison of different parton shower implementations and briefly comment on the 
phenomenology of the colour coherence effect at the Tevatron. 
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Figure 1. The top quark forward-backward asymmetry Afb in 
dependence on the transverse momentum of the pair. Results 
from MC@NLO [1], Pythia [3], Powheg [+], MCfm [: ] and tt+jet 
production at NLO [6] are shown as reported by the D0 and CDF 
Collaborations in Refs. [7] and [S], respectively. 



1 Introduction 

In a number of different measurements the Tevatron Col- 
laborations both CDF and D0 reported results on the top 
quark forward-backward asymmetry that clearly lie be- 
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yond the Standard Model expectations for this quantity, 
see Refs. [7-10]. They did not only determine inclusive 
asymmetries but presented also differential distributions, 
such as ApsifUn), showing the asymmetry in different bins 
of observables O that reflect the kinematic properties of 
the reconstructed top quark system. One such observable 
is the transverse momentum of the pair, for which D0 
made an interesting observation, documented in Ref. [7], 
while studying various Monte Carlo (MC) predictions: the 
leading-order (LO) event generator Pythia [3], using the 
option for approximate angular coherence, displays a qual- 
itatively similar pr.ti dependent asymmetry to that of the 
next-to LO (NLO) matched parton shower, MC@nlo [2]; 
cf. the top panel in Fig. 1. Disabling the coherent shower 
option will restore the expected behaviour for Pythia as 
indicated by the dashed line (magenta) in the figure. CDF 
later on published similar results [ ] confirming D0's ob- 
servation; this is also shown in Fig. 1, at the bottom. 

The definition of ApB used by the Tevatron experimen- 
talists is based on taking the rapidity difference between 
the top and the antitop quark. Ay = y, - yj, and dividing 
the events according to their Ay hemisphere: 



\do)+ \dOl- 



(1) 



The terms (dcr/dO)-t denote the differential (or, if (9 s 1, 
inclusive) cross sections measured for forward/backward 
(+Ay > 0) top quark pair configurations; note that Ay is 
a longitudinally boost invariant quantity. We shall now 
investigate the effects of colour coherent parton showering 
on ApB in more detail. 

2 The QCD colour coherence effect 

We can gain a qualitative understanding by considering the 
large-A^c colour flows in the leading-order partonic qq 
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Figure 2. (Top) Colour flow and QCD radiation in (a) forward 
and (b) backward top quark pair production. (Bottom) Sherpa 
CSshower [11, 12] predictions for all and separately generated 
forward/backward (fwd/bwd) ff configurations at LO. The colour 
coherence effect is illustrated after the first and finally all shower 
emissions have occurred by means of the pr ,,- distribution. 

tt processes, the by far dominant contributions to tT pro- 
duction at the Tevatron. The (leading-)colour flows asso- 
ciated with these partonic interactions stretch from the ini- 
tial to the final state connecting the incoming (anti)quark 
with the outgoing (anti)top quark. They form initial-final 
colour dipoles as shown in the upper part of Fig. 2. The de- 
flection of the top quark out of the incoming quark's direc- 
tion is a measure of the acceleration of the colour charge. 
Colour coherence manifests itself through this deflection 
angle: a mild scatter characterized by a small opening an- 
gle, keeping the forward motion of the top quark, only in- 
duces weak radiation off this qt dipole, see Fig. 2 upper 
left. In contrast, strong scattering and re-direction of the 
top quark leads to a more violent acceleration of colour 
and hence a larger number of potentially also harder QCD 
emissions. Top quark pairs in backward configurations 
therefore experience an increased recoil which we illus- 
trate in the upper right of Fig. 2. Consequently, we find a 
larger number of soft pj,f events correlated with positive 
values of the asymmetry while hard prjr events often con- 
tribute to negative asymmetries - as well known from the 
NLO calculation, cf. Fig. 1. 

We can easily test this rationale in a MC simulation 
which we do by utilizing Sherpa's CSshower implemen- 
tation [11].' We track the evolution of top quark pairs gen- 
erated at leading, fixed order in forward (Ay > 0) and 
backward (A^ < 0) configurations independently. The 



'The use of dipole-like splitting kernels derived from spin-averaged, 
leading-colour reduced Catani-Seymour subtraction terms ensures an ac- 
curate treatment of the soft limit of QCD emission, thus, enables the con- 
struction of a colour coherent radiation pattern. 
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Figure 3. The function F(/3, pj), normalized to its absolute value 
at zero p^; the maximum pj is given as p™™ = ^ 



Figure 4. 

Leading-order QCD 
predictions (calculated 
with MCfm [ ] at NLO 
in tt production) for the 
forward-backward 
asymmetry as a 
function of p'j imposing 
an upper, a lower and 
no bound at all on the tt 
pair mass. 



tilde sign in bdj is used to flag an evaluation at the matrix- 
element level. The results in the lower pait of Fig. 2 clearly 
demonstrate that after the coherent showering phase, ini- 
tial backward configurations yield a harder pj,t spectrum 
than initially forwards moving top quark pairs. The lat- 
ter preferably populate the soft region of pj,^ < 10 GeV, 
and both observations thus agree well with our qualitative 
explanation from above. 

2.1 Analytic approximation of the effect 

The qualitative picture which we argued for in Sec. 2 can 
be made more explicit. Analyzing the real-emission con- 
tribution to the asymmetry, mainly arising from qq — > tig 
processes, one can show, see Ref [ i ], that the differential 
dependence of the asymmetry on the pair transverse mo- 
mentum can be written as 



o-fi dpT 



A|-4 
Nc 



F(J3,Pt) 



(2) 



where pr = pT,tt and p - -yjl - 4mf/s denotes the top 
quark center-of-mass velocity. The Bom and asymmetry 
cross sections are given by &b and (Ta, respectively. 

We can now discuss to what extent a coherent branch- 
ing algorithm reproduces this functional form given in 
Eq. (2). Such shower generators treat the gluon radia- 
tion as coherent emission from the top quarks lines of 
the Born process, therefore overestimate the colour fac- 
tor in Eq. (2) by 2Cf = (A^^ ~ l)/^c (a 60% overesti- 
mate) and approximate the kinematic part of the asym- 
metry amplitude by the coiTesponding Born term times 
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Figure 5. Differential cross sections dcr/dO (in pb/GeV, left panels) and asymmetries Afb(O) (right panels) as a function of the tt pair 
transverse momentum (pr.n in GeV, upper row) and ff pair mass (m„- in GeV, lower row). Various predictions (and their MC errors) 
given by LO event generators are compared to each other. 



dipole-like/eikonal factors. That means we effectively ob- 
tain a description of the asymmetry which is correct in the 
soft gluon limit (apart from the colour factor), but only ap- 
proximate for pt > 0, i.e. in the coherent shower genera- 
tors the coupling- and colour-stripped asymmetry function 
F(J3, pj) given in Eq. (2) can only be "underestimated" by 
F{fi,0) which reads F{fi,0) = -4/3 + Oifi^) when 
Taylor expanded in p. Note that F(J3,pt) becomes less 
negative the further away from the pj - Q limit. This 
is depicted in Fig. 3 for different fixed values of /? where 
we observe larger deviations for increasing top quark pair 
invariant masses. Both the colour factor and kinematic ap- 
proximations lead to a more pronounced AFB(pr,//) depen- 
dence, as seen in Fig. 1 when comparing the Pythia coher- 
ent shower prediction with that given by MC@nlo. 

Finally, Fig. 4 shows the pr > dependence of the 
asymmetry as defined in Eq. (1) which one obtains from a 
full NLO calculation for tf production at the Tevatron.^ It 
is interesting to note that the asymmetry approaches zero 



^The finite part of the virtual correction leads to a positive delta peak 
contribution at pr = which however is not shown here. 



but does not become positive as pr 0. This is because 
the singularity structures present in the numerator and de- 
nominator are different. The denominator diverges faster 
owing to initial-state collinear singularities that cancel in 
the numerator. In contrast to the behaviour at fixed or- 
der, NLO matched or coherent parton showers predict an 
Apsipjii) cross-over, as shown in Fig. 1, at low but non- 
negligible values of pr- This is a consequence of the Su- 
dakov suppression of small pr occurring due to multiple 
soft gluon emission, an effect absent in the fixed-order de- 
scription. This Sudakov suppression yields a spreading of 
the positive asymmetry over a finite region of pj > 0. 

3 Differential asymmetry predictions 

Considering the impact colour coherence has on generat- 
ing differential asymmetries, it is important to study the 
response of standard MC event generators using ff produc- 
tion (without decays) at the Tevatron. For a representative 
collection of parton showers, we present their respective 
asymmetry predictions in Fig. 5 together with the associ- 
ated differential cross sections as functions of pj-,,- (upper 



4 



Monte Carlo event generators & the top quark forward-backward asymmetry 



February 14, 2013 



Table 1. Inclusive Apg (in %) in each shower model, and to 
leading, non-trivial order in QCD (taken from MCfm). 



Model Version Afb Model Version Afb 



row) and m„- (lower row). Far more details regarding this 
MC tools comparison and the choice of parameters can be 
found in Ref . [ ] . Here we only highlight a small fraction 
of the results of this comparison. 

As one can see from Fig. 5, default Herwig-i-+ [13] 
and Sherpa [ ' 2] produce differential Afb sufficiently sim- 
ilar to what one expects from the approximate NLO treat- 
ment mentioned above, especially the rise of Afb with in- 
creasing m,r is remarkable. Using the ot„- observable the 
Sudakov region is stretched/applied over the entire mass 
range which leads to a (an almost entirely) positive Afb 
dependence. Both Herwig-h- (through angular ordering) 
and Sherpa (through dipole showering) have implemented 
QCD coherence in a proper way, however rely on dif- 
ferent strategies regarding initial shower conditions and 
treatment of recoils. This induces differences between 
their predictions as seen in Fig. 5. While the version of 
Pythia 8 [ ] used here does not yet account for QCD co- 
herence effects, hence produces rather small, pj^n as well 
as m„- insensitive asymmetries, the Pythia 6 [3] predictions 
in Fig. 5 are examples for the fact that the older generation 
provides options (tunes) with varying amount of coher- 
ence. The D6T and P(erugia)0 tunes [15] emerged from 
efforts to improve the description of Tevatron and LHC 
data, respectively. In particular the D6T tune amplifies the 
effect of colour coherence as explained in Ref. [ 1 ]. 

4 Inclusive asymmetry and longitudinal recoil 
effects 

Inspecting the numbers in Tab. 1 , it is obvious that colour 
coherent showering not only produces non-trivial differen- 
tial asymmetries but also an inclusive Afb-"* This clearly 
comes as another surprise, but can be explained via event 
migrations arising from longitudinal recoil effects. This 
is specified in Ref. [1] where we also show that the ef- 
fect, already at LO, may give a relevant contribution to the 
inclusive asymmetry. Here, we only want to provide ev- 
idence that migration indeed occurs in coherent showers, 
using the CSshower results displayed in Fig. 6. 

Following the same strategy used to produce the pT,ri 
spectra presented in Fig. 2, it is straightforward to obtain 
the corresponding Ay distributions to examine the migra- 
tion effect directly. A simple visualization of the longitudi- 
nal recoil effect comes in terms of a stretching/widening of 
the gluon emitting qt initial-final dipole. Independently of 
the magnitude of the initial scattering angle, the top quark 
will be slightly pushed forward while the antitop quark re- 
tains its direction such that Ay - H- e is a plausible 

^This can also be anticipated by viewing Herwig++'s and Sherpa's 
positive-valued ApB(m,f) results given in Fig. 5. 
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Figure 6. The Ay distributions for various modes of CSshower- 
ing applied after disjoined as well as joined Ay hemisphere gen- 
eration at leading, fixed order in top quark pair production. 



Table 2. Inclusive ApB results (in %) imposing simple kinematic 
cuts (M = 450 GeV, g = 50 GeV) for the Herwig++ model, the 
CSshower and tt production at NLO (computed with MCi^i). 



Model 


Version 


m,f < M 


> M 


Pt.it < Q 


>Q 


Herwig-h- 


2.5.2 


2.7 


6.0 


5.8 


-14.3 


Sherpa 


1.4.0 


3.5 


9.2 


8.7 


-15.4 


QCD 


LO 


4.1 


9.3 


7.0 


-11.1 



parametrization (e > 0). The results shown in Fig. 6 con- 
firm this conjecture. Several observations are made: (1) 
the imbalance between Ay <Q) and Ay > Q events, i.e. the 
generated asymmetry, is clearly visible by comparing the 
after- with the before-shower result (cf. red vs. black line), 
(2) the migration of initial backward and forward config- 
urations to higher values of Ay produces an overflow of 
events of the former category closing the dip at At/ > 
generated by shifted events of the latter category (cf. green 
vs. blue lines), (3) -H — > - migrations in Ay (i.e. e < 0) 
are suppressed and (4) the largest contribution to the mi- 
gration effect akeady appears with the first emission (cf. 
dashed vs. solid lines). Generally, migrations are small, 
occur locally and clearly favour e > 0. They therefore fuel 
the generation of a positive-valued inclusive asymmetry. 

5 Phenomenological tests 

With our understanding of the origin of the colour coher- 
ence effect, we can study some of its phenomenological 
imphcations which we briefly discuss below. 

5.1 Asymmetry in different phase-space domains 

By focusing on certain phase-space regions, one can am- 
plify the inclusive asymmetry, as shown in Tab. 2. The 
application of cuts isolating large m,j or low pjjt leads to 
an increase of the positive asymmetry. The latter cut is 
particularly useful to efficiently separate off the Sudakov 
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Top quark FB asymmetry vs /! mass @ Tevatron Rmi2 
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Figure 7. ApBO'nn) predictions (center) and ratio plots using colour coherent showering (here computed by Sherpa's CSshower [11]) 
for a prjf domain separation at pj-jj = 50 GeV. A recoil scheme variation (as provided by Sherpa v. 1.4.0) is also shown as well as 
reference curves for lowest-order ttj production employing different scale choices. 
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Figure 8. ApB(pr,,f) distributions obtained after coherent parton 
showering (computed with the CSshower [11]) for the constraint 
of large pair masses, m,f > 450 GeV. The brighter dashed/dotted 
curves indicate the impact of the first two emissions. The green 
solid curve shows the outcome of a recoil scheme variation (as 
provided by Sherpa v. 1 .4.0 [ 1 2]). 



region, i.e. the domain of low px ,, that generates the pos- 
itive contribution to the overall asymmetry. Considering 
the anti-cut (real-emission) region, one can test the most 
striking feature of the transition in pr.n domains: that is the 
sign change in the asymmetry. Because of the approxima- 
tions discussed earlier, coherent showers predict a stronger 
asymmetry flip at large pr^r- For all other cases, they how- 
ever give enhancements similar to the NLO result. 

These single-cut, Afb enhancing effects can also be 
found at the differential level. We exemplify this in 
Figs. 7 and 8 using Sherpa CSshower results. Figure 7 
is shown to compare the different Apsimri) distributions 
associated with the two pj,f domains. It exhibits the sign 
flip in the asymmetry very clearly. The opposite scenario 
is presented in Fig. 8 where small top quark pair masses 
are vetoed and Afb is shown as a function of the pair trans- 
verse momentum. As expected, the increase of Afb is most 
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Figure 9. Apb(j8) distributions obtained from coherent parton 
showering (as implemented in the CSshower [i 1]) for varying 
"top quark" masses, m, and pp collider energies, ^fs. 



noticeable in the first few pj^r bins, those that are close 
to zero. In both figures we also illustrate the impact of a 
recoil scheme variation (as provided by Sherpa v. 1.4.0): 
differences only occur in the high pr^tt region where the 
amount of Tevatron data is too sparse to discriminate be- 
tween the two options. 

5.2 Top quark velocity dependence of the 
asymmetry 

We want to briefly discuss how the coherence effect (gen- 
erating non-zero Afb) changes under variation of the top 
quark mass and collider energy V*- It is thus convenient 
to analyze the /3 dependence of the asymmetry, which be- 
haves roughly like AFB(«^^^)• Figure 9 contains a number 
of Aps(J3) curves obtained from coherent showering using 
different values for the m and s parameters. Apart from 
PDF effects, one finds, as expected, an approximately sta- 
ble asymmetry dependence, if mass and yfs are scaled 
equally. Investigating the limiting cases, it is clear that 
the dramatic increase of Afb occurs as a consequence of 
reaching the collider phase-space boundaries. The contri- 
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bution to the inclusive asymmetry however is small since 
the cross section of the associated phase-space area is tiny. 

5.3 Top quark forward-backward asymmetry from 
the lepton perspective 

We gain information if we define the analysis not only in 
terms of the forward-backward asymmetry, A,j = Afb, but 
also in terms of a lepton-based asymmetry A/ (now includ- 
ing top quark decays). This allows one to exploit the way 
these asymmetries are correlated in the Standard Model 
and beyond; for the details see Ref. [16]. An example 
for coherent shower predictions in terms of (transverse) 
mass observables, see also [17], is given in Fig. 10 relat- 
ing results for reconstructed top quarks (A,,-, A,f (mti), etc.) 
with results where this reconstruction is not (A/, A; (m^^^j^,), 
etc.) or only partially (A„ (//^^), A/(m„-), etc.) required 
(the "O^^" superscript signifies that the transverse observ- 
ables used here are obtained from all final-state visible ob- 
jects and Bj). This introduces additional analysis handles 
and levels, and therefore more directions to cross-check 
the data for consistency or potential biases. 

6 Summary 

Even in the absence of a full NLO treatment, standard 
Monte Carlo event generators can produce significant dif- 
ferential as well as inclusive forward-backward asymme- 
tries in top quark pair production. These asymmetries arise 
from valid physics built into generators relying on (ap- 
proximately) coherent parton or dipole showering. While 
colour coherent shower algorithms cannot be quantita- 
tively correct in every detail, they are able to capture im- 
portant features as known from QCD higher-order calcu- 
lations. Based on these findings, coherent showers are 
qualified to serve as fast tools in estimating the Standard 
Model phenomenology of the forward-backward asym- 
metry."* However, if used in data analyses to determine 
Monte Carlo corrections, one needs to be aware of their 
previously unexpected asymmetry-generating behaviour. 
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■*For a recent example, see Ref. [ ], where coherent shower predic- 
tions are used to estimate the congelation between A„- = Afb and A/, the 
lepton-based asymmetry definable in the semileptonic decay channel. 

^Visit the mcplots . cern. ch webpages for more information. 



Figure 10. Differential forward-backward (A„- = Apg) and 
lepton-based (A;) asymmetries given in terms of various (trans- 
verse) mass observables [ 1 (>]. Here, CSshower results are shown. 
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